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Overview 


*  Background 

*  BHT-600Thruster 

*  Experimental  Work 

*  Collisional  Radiative  Modeling 

*  HPHALL  Simulations 

*  Numerical  /Experimental 
Comparison 
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Background:  Breathing  Mode 


Constant 
neutral  flow 


ion 


0  J^ve 


100 


Seen  through  low  frequency  (10-50k  Hz)  oscillations  in 
discharge  current  (ld) 

Periodic  depletion  &  replenishment  of  neutrals  at  exit1 
Also  referred  to  as  neutral  transit  time  instability-  scales 

with  L channel  neutrals 


•  Previous  time  averaged  measurements  unable  to 
quantify  oscillations  in  plasma  properties 

[1]  Boeuf,  J.  P.;  Garrigues,  L. ,  Journal  of  Applied  Physics ,  vol.84,  no.7, 
pp.3541-3554,  Oct  1998 
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Background:  Diagnostics 


Probe 

•  Measurements 

-  Langmuir-  Te,  density,  e'  EDF 

-  RPA-  ion  EDF 

-  Faraday-  thruster  beam  current 

•  Intrusive-  spatially  limited 

•  Temporally  limited  due  to  sweeps 

-  Lobbia2  (1 0  jis  resolution) 

[2]  Lobbia,  RB  and  Gallimore,  A.D. ,  Rev.  Sci.  Instrum. 
81,073503,  2010 


Emission 

•  Measurements 

-  Line  Intensity  ratio-  Te 

-  Absolute  Intensity  -  Density 

-  Doppler  Shift  -  velocity 

•  Non-intrusive-  capable  of  near  field 
measurements 

•  Line  of  sight  averaging 

•  Measurements  on  ns  timescales 
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BHT-600  Hall  Effect  Thruster 


•  Thruster  tested  w  Xe  at  nominal  conditions 

•  Extensive  previous  experimental  work 

-  Probe-  RPA,  Faraday,  ExB3  6 

-  Optical  measurements-LIF6’7 

[3]  Ekholm  et  al,  JPC,  2006.  [6]  Hargus  et  al,  JPC,  2008. 

[4]  Niemela  et  al,  JPC,  2006.  [7]  Nakles  et  al,  JPC,  2008. 

[5]  Nakles  et  al,  IEPC,  2009. 


BUSEK 


Dimensions 

Rinner  24  mm 

Router  32  mm 
Channel  Depth  10  mm 
Nominal  Conditions 
Anode  Flow  Rate  2.45  mg/s 
Cathode  Flow  Rate  197  pig/s 
Anode  Potential  300  V 
Anode  Current  2.05  A 
Magnetic  Current  2.0  A 
Performance 
Thrust  42  mN 
Specific  Impulse  1650  s 
Anode  Efficiency  55.0% 
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Current  FFT  (A)  Current  (A) 


BHT-600  :  Breathing  Mode  Oscillations 


Time  (us) 


Discharge  Currect  (AC) 
Passive  inductive  probe 
Band  pass:120Hz-  20MHz 


Discharge  Current  FFT 
Spectrum  Analyzer 
with  FFT  averaging 
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Timing  System 
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Relative  Signal 


Timing  System 


1 — | — I — I — I — I — | — TIT — I — I — | — I — I — I — I — | — I — I — I — I — | — I — I — I — I — [ 


Gate  Delay 

(13(18) 
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Gate  Width  (Ijlis) 


i  i — r 

Current 


Triggering  Circuit 
Camera  Gating 


Trig.  Circuit  Component  Delay  (2.2p,s) 


Time  (us)  Shots  accumulated  over  exposure 

(0.5s),  ~19K  shots  /measurement 
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System  Schematic 


'Chamber  1 ,  AFRL  Edwards  AFB 

Stainless  Steel 
Pumping  Speeds 
of  48,500  L/s  o 

P=  9. 1  (jTorr  AA?™ 


Thruster  Current 


Thruster 


Collection 
Volume 

6  mm  diameter 
fiber  coupler 


200jim 

fiber 


1 .25  m  Spectrometer 
300  G/mm,  500  nm  blaze 
Dispersion=  2.6  nm/mm 


"V" 


Passive  Inductive 
Probe 


Voltage  Signal 
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Collisional  Radiative  Modeling  (CRM) 


*  Predicts  emission  by  modeling 
collisional  excitation  and 
allowed  radiative  decay  paths 

*  KCD8  Metastable  Modeling 

—  Treated  as  virtual  ground 
—  Assumed  in  equilibrium 

Simplified  Xe  Collisional 

Excitation  Processes 

e~  +  Xe  -»  Xe*  +  e~ 

Xe+  +Xe  ^  Xe*  +  Xe+* 

Xe2+  +  Xe  ^  Xe*  +  Xe2+* 
e  +  Xe  m  Xe *  +  e 


Simplified  KCD  Diagram 


ground 


[8]  Karabadzhak  et  al.,  Journal  of  Applied  Physics ,  2006 


Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


li 


CRM:  KCD  Model 


E  x  B  probe 
measurements 


*Metastable  rate 
approximations 


he 

hel(X)=—{N()Ne) 


kx 


eO 


+  aki0  + 


1  —  a 


kA 

k20 


+ 


/ 

f(Te,  a)^0.01%  —  0.3% 
^Equilibrium  assumption 


•  Ions-  uniform  velocity,  LIF  Empirical  excitation 

•  e-  -  Maxwellian  EDF=f(Te)  cross  sections910 


hel  (A  l) 

hel  (^2) 


fia.u^Tg) 


[9]Chiu  et  al,  Journal  of  Applied  Physics ,  2006.  [1 0]  Sommerville  et  al,  Journal  of  Prop.  &  Power,  2008. 
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1(823/828) 


CRM:  KCD  Model 
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Normalized  Emission  Measurements 


0  5  10  15  20  25 

Time  (us) 


•  Intensity  in  phase  with  discharge  current  oscillations 

•  Similar  visible  emission  fluctuations  seen  by  Liu  et  al.11  in  thruster  high  speed  imaging 

•  Small  phase  shift  (~2  jus)  seen  between  intensity  and  discharge  current 

-  Corresponding  to  a  8  km/s  electron  axial  velocity 

-  In  agreement  with  5-10  km/s  electron  axial  velocity  predicted  by  HPHall12 


[11]  Liu  et  al.,  IEEE  T.  Plasma.  Sci.,  (in  press)  [12]  Scharfe,  M  K,  Koo,  J  W,  personal  communication,  2011 
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Electron  Temperature 
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•  Line  ratio  method  divergence  in 
cooling  phase  of  the  cycle  may  be 
due  to: 

-  Non-maxwellian  EEDF 


•  Low  emission  signal  (SNR=1 .9  dB) 
in  low  ld  portion  of  cycle 

-  Higher  uncertainties  in  Te 
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Electron  Temperature 


^823/^828 

^834/^828 

Te 

eV 

6.6  ±  0.6 

6.2  +  1.8 

Te/Te 

% 

27  ±4 

38  +  15 

ps 

-4 

-6 

Comments 

•  High  SNR 

•  Lower  Te  uncertainty 

•  Low  SNR 

•  Higher  Te  uncertainty 

•  Independent  of 
metastable  approximation 
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HPHall:  Overview 


*  Radial-axial  hybrid  particle-in-cell  (PIC) 

—  Fluid  electrons 
—  PIC  ions  and  neutrals 
—  Quasineutral 

—  Electron  mobility-  axially  varying  effective 
mobility 

—  Simulations 

—  Time  step  =  0.25  ps 
—  Varying  Inverse  Hall  Parameter  KB/1 6 
—  Properties  evaluated  at  4  axial  locations: 

A)  Channel  Near  Anode  C)  Near  Plume 

B)  Channel  Near  exit  D)  Far  Plume 

—  Experimental/numerical  comparison: 

—  Time-averaged  quantities 
—  Breathing  mode  behavior 
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HPHall:  Case  Comparison 


Units 

Measured 

HPHall 

Case  1 

HPHall 

Case  2 

HPHall 

Case  3 

HPHall 

Case  4 

T 

mN 

39 

38 

38 

38 

41 

Isp 

s 

1530 

1530 

1510 

1530 

1660 

f  BM 

kHz 

38 

34  +  5 

54  +  5 

31+5 

22  +  5 

Id 

A 

2.05 

2.18 

2.21 

2.20 

3.64 

Id/Id 

— 

32% 

7  % 

11  % 

7  % 

6  % 

Te 

eV 

6.6  ±0.6 

12.4 

10.4 

12.7 

25.8 

Te/Te 

— 

32  ±  8  % 

4% 

6% 

4% 

2% 

•  Accurate  thruster  performance — within  10% 

•  Difficulty  with  predicting  breathing  mode  behavior 

-  Oscillations  in  ld  and  Te  significantly  lower  than  observed 

-  Higher  than  observed  Te-  possible  mechanism  to  increase  mobility 

•  Inverse  Hall  parameter-strong  influence  on  breathing  mode 

-  Variation  in  breathing  mode  frequency  and  amplitude  of  oscillations 
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•  Accurate  thruster  performance — within  10% 

•  Difficulty  with  predicting  breathing  mode  behavior 

-  Oscillations  in  ld  and  Te  significantly  lower  than  observed 

-  Higher  than  observed  Te-  possible  mechanism  to  increase  mobility 

•  Inverse  Hall  parameter-strong  influence  on  breathing  mode 

-  Variation  in  breathing  mode  frequency  and  amplitude  of  oscillations 

•  Further  analysis  based  on  Case  1 
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HPHall  Results:  Time-Resolved 


•  Te  -  Slight  correlation  with  ld  seen 
in  plume,  and  near  anode 

-  Out  of  phase  with  discharge  current- 
similar  to  experimental  observations 

•  N0  -  Little  correlation  with  ld 
except  for  near  anode 


•  Nj/Ne  -  -  strong  correlation  with  ld 
at  all  axial  locations 

-  Oscillations  nearly  in  phase  with 
discharge  current  oscillations 
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HPHall  Results:  Cross  Correlation 


Te  -  Slight  correlation  with  ld  seen 
in  plume,  and  near  anode 

-  Out  of  phase  with  discharge  current- 
similar  to  experimental  observations 

N0  -  Little  correlation  with  ld 
except  for  near  anode 


Nj/Ne  -  -  strong  correlation  with  ld 
at  all  axial  locations 

-  Oscillations  nearly  in  phase  with 
discharge  current  oscillations 


Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


21 


HPHall  Results 


500  550  600  650  700 
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HPHall  Results 


Axial 

Location 

X 

x/x 

freq 

(kHz) 

^"delay 

(US) 

A 

T1d,X 

(degrees) 

R 

Discharge  Current 
\ 

2.2  A 

7% 

34 

- 

- 

1.00 

A 

9.5E+18  m  3 

4% 

33 

6.4 

77 

0.75 

Neutral 

Density 

no 

B 

2.1E+18  m  3 

5% 

33 

-22.9 

-275 

0.26 

c 

5.6E+17 m  3 

7% 

28 

7.1 

73 

0.11 

D 

3.9E+17  m  3 

7% 

33 

5.6 

68 

0.09 

A 

1.2E+24  m  3 

8% 

33 

1.9 

23 

0.90 

Ionization 

B 

1.7E+23  m  3 

9% 

34 

1.1 

14 

0.75 

Rate 

n, 

C 

1.9E+22  m  3 

12% 

34 

0.4 

5 

0.76 

D 

6.3E+21  m  3 

13% 

33 

0.4 

5 

0.76 

A 

1.1E+18  nv3 

7% 

34 

0.4 

5 

0.91 

Electron 

B 

6.4E+17  nv3 

8% 

34 

0.4 

5 

0.93 

Density 

ne 

C 

3.5E+17  nv3 

8% 

34 

-1.1 

-14 

0.90 

D 

2.4E+17  nv3 

9% 

33 

-1.1 

-14 

0.89 

A 

27.3  eV 

2% 

34 

-7.1 

-87 

0.60 

Electron 

B 

23.0  eV 

2% 

34 

-11.6 

-143 

0.31 

Temperature 

T 

C 

16.5  eV 

4% 

34 

5.6 

68 

0.48 

e 

D 

12.5  eV 

5% 

34 

4.9 

59 

0.58 
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Summary 


•  Time  resolved  emission  measurements: 


•  Emission  of  823,  828,  and  834  nm  lines 
found  to  fluctuate  nearly  in  phase  with  ld 


•  Te  fluctuations  of  32%  ±  8%,  -5  ±  1  ps  out  of 
phase  with  ld 

Units 

Measured 

HPHall 

Case  1 

•  HPHall  simulations  run  with  varying 

T 

mN 

39 

38 

effective  mobility  profiles: 

•  Profile  choice  strong  influence  on  breathing 
mode  oscillations  (frequency,  amplitude) 

•  Little  effect  on  overall  performance,  all  cases 
w/in  1 0%  of  observed  values 

hp 

s 

1530 

1530 

f BM 

kHz 

38 

34  +  5 

h 

A 

2.05 

2.18 

Id/Id 

— 

32% 

7  % 

Te 

eV 

6.6  +  0.6 

12.4 

Te/Te 

— 

32  +  8  % 

4% 

•  HPHall  Breathing  mode  predictions: 

QlJe 

deg 

68°+  14° 

68T 

•  Significantly  under  predict  oscillation 
magnitude  for  both  ld  and  Te 


•  Higher  than  observed  Te  values  maybe  be 
result  of  need  for  increased  mobility 

•  Phase  shifts  found  to  be  consistent  with 
current  observations 
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